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A Lindblad master equation for a harmonic oscillator, which describes the dynamics 
of an open system, is formally solved. The solution yields the spectral resolution of the 
Liouvillian, that is, all eigenvalues and eigenprojections are obtained. This spectral 
resolution is discussed in depth in the context of the biorthogonal system and the 
rigged Hilbert space, and the contribution of each eigenprojection to expectation 
values of physical quantities is revealed. We also construct the ladder operators of 
the Liouvillian, which clarify the structure of the spectral resolution. 
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I. INTRODUCTION 



The time evolution of an open quantum system interacting with external quantum sys- 
tems is non-unitary. Therefore we cannot apply directly the usual quantum mechanics to 
the open systems. Instead, we think the total system which is the joint system of a tar- 
get system and an environment to be a closed system to which the quantum mechanics is 
applicable, and extract the time evolution of the target system by taking the partial trace 
over the environment. It is called the Reduced dynamics. Kraus^ showed that if there is 
no initial correlation the time evolution of a reduced density matrix is complete positive. 
Lindbladii and Gorini et al.- obtained the general form of the generator of the complete 
positive djTiamical semigroup (Liouvillian) £, which describes the time evolution of density 
matrices. A master equation with this generator is called a Lindblad master equation. 

In this work, a system of a harmonic oscillator is considered, which interact linearly with 
a heat bath of bosonic field, an environment. The Markov approximation and the van Hove 
limit yield the Lindblad master equation^>i^ii^ii^ of the form 

p = Cp = —i[H, p] + 7 ^apa^ — -{a'^a, p}^ + 7' ^a^pa — -{aa\ p}^ 

H = ooa^a, - = e"^". 
7 

Here, let Jif be the separable Hilbert space of the states in the target system. H is the 
Hamiltonian of the system, u G ]R>o a angular frequency, and a\a are respectively the 
creation and annihilation operators, which are densely defined unbounded operators on Jif 
and satisfy the commutation relation [a, a''"] = 1. We can choose a set of the eigenstates of 
the Hamiltonian { I n)}„g2>o) H\n)=ijjn\n), \n) = \/n + 1 |n + 1) , a\n) = ^Jn\n — \) 
as a complete orthonormal system (C.O.N.S.) of 7, 7' G ]R>o are related to the damping 
rate and /3 G M>o is an inverse temperature. A density matrix (state) is a trace class operator 
on i.e. p G Ti(^), which satisfies Trp = 1, (?/^| p \^) > 0, V G ,^ . Therefore ^ is 
a differential equation on Ti(t^). But we will restrict the space of density matrices later, 
because we treat the unbounded operators as observables. 

This equation is well-known, and its derivation can be found in, e.g., Ref.-. As an 
application of this equation, for example, the analysis of the cavity quantum electrodynamics 
may be considered. The goal of this work is to solve this equation in a way to clarify the 
structure of the Liouvillian L. 
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There have been trials to solve this equation (see Refs.-~-ii^). Briegel and Englert have 
obtained the eigenvalues and the biorthogonal eigenbasis in terms of creation and annihila- 
tion operators by using normal ordering^, and Tay and Petrosky also have done in terms of 
the Wigner function^^. But in this work we derive the solution by eigenprojection method 
in terms of creation and annihilation operators without normal ordering, and present the 
spectral resolution of the Liouvillian C If £ has an spectral resolution 

where Am is m-th eigenvalue and Ilm, is corresponding eigenprojection, the solution of the 
master equation is 

p(t) = e*V0) = 5^e^'"*nmp(0). (3) 

m 

The paper is organized as follows. In section [ITl we derive the solution of ([T]) in the form 
of ([3]), and obtain the presumed eigenvalues and eigenprojections. In section UTTl we confirm 
that they satisfy the property of the eigenprojections of the Liouvillian ([2]). These eigen- 
projections have the structure of the biorthogonal system. Although the right eigenvectors 
are in the left eigenvectors are unbounded operators. This is due to the unbounded 

Hamiltonian. In section IIVI we give a mathematical justification of this spectral resolution 
on the basis of the rigged Hilbert space (Gel'fand triplet). In section |V] we discuss about 
the time behavior of this solution. In section IVII we calculate the time behavior of the ex- 
pectation values and variances of some physical quantities exactly. Although the number of 
eigenprojections is infinite, we need only finite number of eigenprojections. Finally, in sec- 
tion IVIII we show the ladder operators of the Liouvillian. The eigenvalues of the Liouvillian 
have two indices, and there are two sets of ladder operators. The mathematical notions used 
in this paper are summarized in Appendices. 



II. DERIVATION OF A FORMAL SOLUTION AND 
EIGENPROJECTIONS 

In this section, we acquire the general solution of ([1]) formally. First to eliminate the 
commutator and anticommutators, we do the transformation 

pit) = ^i-i^N-,'/2)tpl^yi,n'N-,'/2)t^ ^ ^ _ ^ ^ I±l^ N = a^a. 
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Then ([T]) becomes 

P = ^e-^^'apa) + ^e^^'a^pa. (4) 
To solve this equation we introduce the superoperators /Ci, /C2, /C3 whose actions are 

ICip = apa\ IC2P = pa, JCsp = {N + 1/2, p} . (5) 
These superoperators have the algebraic structure of s[(2). 

[/Ci, /C2] = /C3, [/C3, /Ci] = -2/Ci, [/C3, /C2] = 2/C2 
Existence of this algebraic structure may allow us to expect to find a solution of the form 

In terms of the superoperators, the time derivative of p becomes 

P = {/ie-2/3/Ci + (/2 - 2/2/3 + /le-^^Vl) ^2 + (A - /ie-^^/2) } p (6) 

Here we use the Hadamard lemma, i.e. for the elements X, Y of Lie algebra g, 

e^Xe-^ = X + [F, X] + i [Y, [Y, X]] + • • • = e^^^^^X, 

where ad(y)X = [V, X]. On the other hand, we can write down the right hand side of (j4]) 
in terms of superoperators 

p = 7e-2^*/Cip + 7'e'^*/C2P (7) 

Comparing ([6]) with (^^, we notice that fj{t) (j = 1, 2, 3) satisfy the simultaneous differential 
equation 

/3 - /ie-'^V2 = 0. 

We solve this equation with the initial conditions /j(0) = (j = 1, 2, 3). First we can extract 
the differential equation of /2 



This is a Riccati equation. Once we obtain the solution of /2, we can solve the other 

equations just by integration. The solutions are 

1 _ p~h-i')t 

f^it) = : ^ -, (8a) 

7 

V 1 _ p-(7-7')* 

7 

/e(7-7')t _ V\ 

Consequently, we obtain the formal solution of ([1]) in terms of superoperators 

p{t) = e(-*f^^-V/2)t|g/2(t)/C2g/3(i)/C3g/i(i)/Ci^(^Q^)|g(in*7V-7'/2)t^ 

By using the following relations 

^(_injV-7'/2)t|g/2{t)/C2^|g(in*7V-772)t ^ g-7'ig/2(t)e-2^'/C2|g-if^iVt^gin*7Vt| 
pfzmzp ^ g/3(t)(7V+l/2)^g/3W(A^+l/2)^ 

([5]) is rewritten as 

(10) 



p(t) = /(x)e-^"^(")^^e^^^ <^ (e-'"*a;^/(x)) (e-"^(")'^ie'^V(0)) (e'"*a;^/(a;)y 



1 - 21 

x(t) = e-(^-^')*, fix) - ^ 
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Here, we describe as a matrix in the basis of eigenstates and introduce some formulae 

oo 

= n \n) 



(n\ 



n=0 



|0) (0| {e^^p} |0) (0| = |0) J2 (^1 P 1^) (0| = |0) (0| Trp, 

n=0 

|0) (0| = ^e"'"^ |n) {n\ = e'^^'''^ 

n=0 

The superoperator e^^ brings about essentially the trace operation, and acting e'^^ on |0) (0| 
makes the equilibrium density matrix. Using above formulae. 



p{t) = Y J_ -i{m-n)ut^n^ j,^^^m+n+l^-^xmK, 

^-^ mini 

m,n=0 

X ((at)'"e-^^'^^''a") Tr {a™ (e-^^(^)'^V(O)) (a^)"} 



mlnlplql \ 7 

m,n,p,q=0 ^ 



x(a^)™+Pe~^""^"a"+PTr {a"^+V(0)(a^)"+''} • 



Expanding f{x) with the formula 



1 X 
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-s \ ^ / I \ r 



^ r!s-l! Vt 



we obtain the formal solution 



m,n,p,q,r=0 



{■m + n + p + q)\m\n\p\q\r\ 



,\m+n+p+q+l / i\q+r 

1 _ 21 \ ' ' 
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(11) 



X e-('^+P+5+^0(^-^')*e-^(™-")'^*(a^)'™+Pe-^""^'^a"+fTr {a"^+V(0)(a^)"+'^} . 
Here we notice that the time dependent parts in (fTTl) are only 



exp 



m + n 



+ p + q + r\ (7 — 7') — " r ^ 



Assembling the combinations (m, ?7,,p, g, r) which have the same exponent, we can write the 
solution in the form of (|3]). For example, 



0,0,0,0,0 



1,0,0,0,0 



0,1,0,0,0 



2,0,0,0,0 



0,2,0,0,0 



1,1,0,0,0 



0,0, 1,0,0 



0,0,0,1,0 



0,0,0,0,1 



Ao,o — 0, 



Ai,i = 
Ai,_i 



7-7 



7-7 



lUJ, 



no,op 
ni,ip 



1 - e-^^^^'^Trp, 
7/ 



ni,_ip = ( 1 - 21 ) e-''""''^aTr {pa^} 
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1 y 



A2,2 = -(7 - 7') - 22u;, n2,2P = - (^1 - ^ ) (at)2e-^-«^'^Tr {aV} 



A2 -2 = -(7 - 7') + 2iw, Ha _2P 
A2,o = -(7 - 7')) n2,oP = 



7 \ „~/3a;ata^2 



7 



Tr{p(at)2} 



1 1 e ""a" 

^] a^e-^'^^^^'aTr {apat} 

llate-'^-'^^'^aTrp 
7/ 7 



1--) e-'^'^'^^'^Tr {apa^} 



+ 1 



7 

H ) He-^'^'^^'^Trp. 
7/ 7 

Executing this program, we presume that the eigenvalues are characterized as 



Ai,fc = "2(7 - 7') - ^ku, j e Z>o, A; = -J, -j + 2, ■ ■ ■ , j - 2, j. 



(12) 
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which agree with those in Ref.— , and the corresponding eigenprojections are, for k > 



(13a) 



^'''^ ~ {k + m)\{k + n)\{l - m)!(Z - n)\m\n\ \ ~) \~ 

X (a^)'=+"^e-^'^'^^"a"^Tr {a'^+Xa'^)"} , 

^ (-ir+"(fc + /)!/! / _ yy+n.+n-M .^W^n 

^^^{k + m)\{k + n)\{l - m)\{l - n)\m\n\ \ J V7/ 

]\m ~l3uja''a k+mrY ( nf n^'\''+^ 



X a' e ^ a 



Tr{aXa1')^'+"} . 



Eventually we have obtained all the presumed eigenvalues and eigenprojections. 

III. PROPERTY OF THE EIGENPROJECTIONS 

In this section we will show that the eigenvalues and eigenprojections which we obtained 
in the last section certainly satisfy 



Hj- fcllj/ fc/ = 5j fc/IIj (orthogonality and idempotence), 
Ilj = 1 (completeness). 

First, we notice that the eigenprojections can be divided into the operator part and 
the trace part. We expect that the eigenprojections have the structure of biorthogonal 
system. To explain this idea, we introduce the bra-ket notation for the vector space of 
operators. We describe a Hilbert-Schmidt operator A G T2(^) as a ket \A)). Since %2{J^) 
is a Hilbert space, we can also treat B"^ G %2{'^) as a bra {{B\. We define the Hilbert- 
Schmidt inner product between bras and kets {{B\A)) = TtB'^A. When we restrict the 
ket space to %i{J^), the bra space becomes 'tB{J^). The pairing is defined similarly as 
{{B\A)) = TtB^A, a G Ti(^), B G <B(^). Let O be an superoperator on the ket 
space. The adjoint superoperator (9^ is the operator defined on the bra space whose action 
is defined as, {{O^B\A)) = {{B\ O \A)). 
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Now we get back to the discussion of the eigenprojection. For > 0, we define 



I 



\j,k)) ^ v^5^ Afc,,,„(at)'=+"e-'^-'^''^a-, / = LJi, (14a) 



m=0 
I 



{{j, k\ ^ v^$^ i?;.,z,„(a^)"a'=+", (14b) 



n=0 



Ak,l,m= ^ \,./ 1--) , (14c) 

[k + m)!(/ — m)!m! V 7 / 

B.,.^ ^-!r d-^rm'". 



(A; + ri)!(/ — n)!?!! \ 7/ \7 



,N fc+l / /N / 

7\ /7^ 



C,,^{k + l)m^l--j ^-j. (14e) 

I J, fc)) is shown to be the right eigenvector of the Liouvilhan corresponding to the eigenvalue 
Xj^k, and {{j, k\ to be the left (generalized) eigenvector. Then the presumed eigenprojection 
can be written as 

U,,k=\j,k)){{j,k\. (15) 

For k < these eigenvectors are defined similarly. Consequently we obtain the spectral 
resolution of the Liouvillian as 

oo j 

^ = b'^' - 20) - 2/ 1 . (16) 

j=0 1=0 

It is easy to see that (a^)'^e~^'^"^"a'^, 171,11 & Z>o is a bounded and trace class operator, 
whose trace is 

Tr {a-(at)"e-^-«^«} = 6^^^^^-^^. (17) 

This formula is derived easily by induction. Therefore the right eigenvectors belong to the 
space of trace class operators %i{J^). This is desirable because axiomatically the space of 
the density matrix is in %i{J^). However the left eigenvector is obviously unbounded and 
not in %i{Jif). This problem is treated in section HVl 
Now what we should prove is 

^ \j, k)) = Aj- fc \j, k)) , ((j, k\C = \j,k {{j, k\ , (18a) 
{{j,k\j',k')) = 6j,,y6k,k' (orthonormality), (18b) 
|j, k)) ((j, k\ = 1 (completeness). (18c) 
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Let £^ be the adjoint operator of the LiouviUian with respect to the Hilbert-Schmidt 
inner product, 



= i[H, p] + 7 ^a^pa — -{a^a, p}^ + 7' (^pa) — -{aa\ p} 



For k,m,n & Tj 



>o 



£{(at)'=+'"e-^'^'^^"a™} = !^-ikuj - ^{k + 2m)(7 - 7')| (a^)'^+'"e-'^"^"^"a™ 
£{(a^)™e-^""^"a^+"} = - ^{k + 2m)(7 - 7')| (at)"^e-^""^"a^+'" 

£t{(at)^+"a"} = jzfcu; - ^{k + 2n){-f - 7')| (a^)'=+"a" + 7'(fc + n)n(a^)^+"-^a"-^ 



Using above results, we can show f llSal) . 

Next we confirm the orthonormahty. Althogh the orthogonahty is trivial from the prop- 
erty of the eigenvectors, we check this by explicit calculation. 

I I' 

{{f,k'\j,k)) = ^Ck,iCk',i' 5^ 5^ Afc,,„i?fc,,,,„Tr |a™(at) V+"(at)^+-e-^'^"''^} (19) 

m=0 n=0 

Apparently the trace vanishes unless k = k'. For h G Z>o, 



a"^(af)"a^+"(at)'=+"^e-'^'^'^^'^ \h) = {h + m - {n - l)}{h + m - {n - 2)} ■ ■ ■ {h + m) 

x{h + l){h + 2)---{h + k + m)e~^^'^ 



Setting ^ = e 
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Using above equation, 

dC 

dC, dC, dC, 



dU \dL 



/ J \ " £•? 



dU (1-0' 

I I An min(m,n) , \ i / -i ^ 

m!n!4 x ^ [k + m + n — ajl f 1 — E, 

n+l 2-^ 



Substituting this trace to (flQll 

' (1 - 0^=+^ £^^{^^{k + my.{k + n)\{l- m)\{V - n)! 

min(m,n) , \i /i ^ 

(fc + m + n — aj! / 1 — 4^ 



(-l)'"(fc + m + n-a)! 



n=0 

; min(m,n) 



^-^ ik + m)\{l — m^lim — a)\{n — a)\a\ 

m=0 a=0 ^ ' ^ ' ^ ' ^ ' 

If we analyze the above sum for / > it suffices. We claim that 

^-^ ^-^ [k ^ m]\[l — m)\[m — a]\[n — a]\a\ II 

This immediately leads to the orthonormality (118bp . 

Now we prove (l20l) . Expanding (1 — and changing the order of the summation, the 
left hand side of fl^U]) becomes 

' "^^^^'"^ A (-ir+"+^(A; + m + n-«)! , 

h h (k + m)\il - m)\im - a)\in - a)\ia - 

_ V(-1)'.^-/3V (-ir V i-ink + m + n-a)\ 
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Here we introduce a formula. For p,q,r & Z>0) P ^ Q ^ f, 

E(-i)^r-^^^4ir^^ = (-i)'^-O' (21) 

^ [p + s - qy.{q - sy.sl 

which is proved in Appendix We apply this formula to the summation about m by setting 
p = k + 1, q = I — a, r = I — n, s = m — a. Futhermore we introduce another formula. For 

p,r e Z>o, p >r 



Applying this formula to the summation about a, we immediately obtain the right hand 
side of ( !20l) and finish the proof of the orthonormality. The formula ( !22|) is heavily used in 
section IVll 

Finally we prove the completeness. If we ignore the left eigenvector problem temporarily, 
it is sufficient to prove that Span ({|j, k))}j^k) is dense in Ti(J^) and use the logic in appendix 
B. We see that 

Span({|j, A:))},,,) = Span ({(at)-e-^-»^'^a"}-„=o) = © 

by the recombination of the basis. We will construct {|m) (^|}mri=o; "which is the C.O.N.S 
of the Hilbert-Schmidt space T2(^), as the limit of the point sequences on &. We use the 
next formulae^ 



■.e-^''':= |0) (0|, 

where : A : represents the normal ordering of the operator A. Setting r = e'^'^ — 1, 



oo p 



e =: e := > — a' ra 



p=0 ^ 



oo 



-1)P 



: e"'^'" := ^ ^-^{a^)PaP = |0) (0| . 
p=0 ^' 

By acting a\a on the above equations, 

oo p 

(at)''+"^e-'^""'"a'?+" = (a^)'? ^ ^(a1")*'+™a^'+"a 

p=0 ^' 

\m) {n\ = ^=y^-L(aY 
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p+m^p+n 



Then we can obtain \m) (n\ from {(a^)™'e ^'^"^°a"}J^„^Q as the following way, 



ym\n\ r. Q.- ^ P- 

* g=0 p=l 



We have constructed the C.O.N.S. of T2(^) as the limit of the point sequences on &. 
Therefore Span ({|j, is dense in %2{'^) and dense in Ti(J^). The proof of the 

completeness taking the left eigenvector problem into account will be done in section IIV[ 



IV. RIGGED HILBERT SPACE OF THE OPERATORS 

In this section we discuss the space which the eigenvectors of the Liouvillian belong to 
more rigorously. The right eigenvectors are trace class operators, but the left eigenvectors 
are unbounded operators and not in %i{Jf). To treat this problem we need the concept of 
the rigged Hilbert space. The Hilbert space which we deal in is the Hilbert-Schmidt space 
'2^2 (^)- 2^1 (^) is a Banach space with respect to the trace norm and dense subset of 
%2{^)- Therefore 

Ti(^) C M-^) C Ti(^)' = *B(^) 

form a rigged Hilbert space. However the dual space is the space of bounded operators. 
We need to think about more restricted space 

= {pe l2{^)Kp{a^T e Ti(jr), Vm,n G Z>o}. (23) 

with the seminorms 

imU,„= l|a'"A(atr||i, m,neZ>o. 

This space becomes a Frechet space. For example (a^)™e~^'^'*^"a'^ and \m) {n\ belong to 
D{J^). The right eigenvectors are in includes the C.O.N.S of T2(^), thus 

D{J^) is dense in %2{'^) with respect to the topology of %2{'^)- If we think about the 
triplet 

c T2(=^) C (24) 

D{M')' includes the operators like (a^)™'a", which appear in the left eigenvectors. The left 
eigenvectors are in D{M')' . 

Now we consider the spectral resolution of the Liouvillian using the logic in appendix [Bl 
{|j, k))}j^k is the set of the eigenvectors of the Liouvillian and a complete system of S)(J^). 
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From (llSbp . {\j,k)) , {{j,k\}j^k is a biorthogonal system. Therefore {{{j,k\}j^k is the set of 
generahzed eigenvectors and complete in generahzed meaning. Consequently the spectral 
resolution ( TTB]) is justified. 

Axiomatically, a density matrix is a normalized, positive semi-definite trace class operator, 
but above discussion imply that we should restrict the class of density matrix more strictly 
if we treat the unbounded operators as observables. For example, 

n=l 

is a density matrix, but the energy expectation value of this density matrix is divergent. 
Therefore this restriction is physically reasonable. 



V. TIME BEHAVIOR OF THE SOLUTION 

Eventually we have obtained the solution of the Lindblad master equation for an harmonic 
oscillator ([1]). 

oo j oo j 

j=0 1=0 j=0 1=0 

This solution shows the time evolution of each eigenprojection explicitly. e~2(')'~'>'')* has the 
effect of dissipation and e"*^'^*, k = j — 21 represents the phase rotation, whose role become 
clear in section IVII Thus the eigenvalues and the eigenprojections are classified according 
to the dissipation quantum number j and the phase rotation quantum number k. The only 
eigenprojection which belongs to j = and does not decay is 

no,op(0) = (^1-^^ e-^^"^"Trp(0) = (l - e"^^) e-^'^^^". (26) 

This is nothing but the density matrix of the equilibrium state peq- Because other eigenpro- 
jections decay, the state converges to the equilibrium state as t — )■ oo. 

Next we confirm that this solution is appropriate for the time evolution of a density 
matrix. First we can show that only IIo^oP have nonzero traces. Apparently only A; = 
eigenprojections has nonzero trace, therefore we study the case j = 2/, A; = 0, / G Z>o. 

((lb,0)) = VC^^i:M,^Jr{iaYe-^-^'^a-} = -^^^ti'^)^ = ^'^o- 

m=0 <! I 1 — ^ I m=0 
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This leads to the trace preserving property Trp(t) = Trp(O) (Vt G M>o), and the conservation 
of probabihty Trp = 0. Self-adjointness of this solution is obvious. It is difficult to prove the 
positive-semidefiniteness by the explicit calculation. However the Liouvillian is a generator of 
complete positive dynamical semigroup, and the solution of Lindblad equation must preserve 
the positivity. 



VI. CALCULATION OF SOME PHYSICAL QUANTITIES 

In this section we calculate the expectation values and the variances of some physical 
quantities. The formula f l22p . which is introduced in section UTTt plays an important role 
in the following calculation. First we calculate the expectation value of energy {E)^ = 
TT{pit)H}. 

oo j 

{E), = a;^^e-^(^-^')*e-(^--2')-* ((ata|j,j - 2/)) - 2/|p(0))) 

j=0 1=0 



J]e-'(^-^')*((ata|2/,0))((2/,0|p(0))). 



1=0 



Here we calculate ((a'''a|2/, 0)). 



((ata|2/,0)) = v^$^^o,z,mTr[, 



m=0 



7 



-2 I 



m=0 



^! , 7' 

Ti \V~\^^ + ~ 

[I — m\)ml 7 



-2 



7 



Thus the expectation value of energy is 



1_ 7 

7 



7 - e 



-(7-7')t 



1-1 

7 



_7 

7 



^/,0 

7 



7 



UJ 



6i,o--ki ((2/,0|p(0))) 

7 



Tr {p{0)a^a} 



(27) 



where {■)^^ means the expectation value at the equilibrium state, that is Tr {pe^if}. The 
convergence {E)^ — > (E)^^ , t — )■ oo is consistent with the convergence of the state. This 
result agrees with that in Ref.-. We notice that because of the formula (|22|) only the 
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eigenprojections {j, k) — (0, 0), (2, 0) contribute to the calculation. Similarly the variance of 
energy is 



In this calculation the eigenprojections (j, A;) = (0,0), (2,0), (4,0) contribute to the second 
moment of energy. In general the eigenprojections {j,k) — (0, 0), (2, 0), • • • , (2p, 0) con- 
tribute to the p-th moment of energy. As expected, these results describe the convergence 
from the initial values to the equilibrium values. If we know the initial and equilibrium 
energy expectation value, we can calculate the time behavior of energy expectation value 
exactly. 

Next we calculate the expectation values and variances of position and momentum. The 
position and momentum operators are described in terms of creation and annihilation oper- 
ators. 



+ e-^(--^')* (2 {E)l - 4 {E)^^ {E), - u {E), + {E\) 
{AE,r = {E\ - {E)] 



(28) 




(29) 



X — 




a — a 



The expectation value of position is 



{x) 




1 



t 



1 



1 




V2u 



e-^(7-7')« [e-^'^*Tr {ap(0)} + e^'^^Tr {aV(0)}] 




(30) 



and that of momentum is 




v^e-5(T-7')tini [e-^'^*Tr{ap(0)}] . 



(31) 
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These results also describe the convergence from the initial values to the equilibrium values. 
The {m,n) — (1, 1), (1, —1) eigenprojections contribute to the expectation values. Similarly 

- i^X + [^"''"Tr {aV(0)} + e^^-Tr {(at)V(O)} 

+Tr {a^ap(O) } + IV {aaV(O) } - 2u {x^)^^ 
^\ - - [e---Tr {aV(0)} + e^^Tr {(at)V(O)} 

-IV {atap(O)} - Tr {aaV(O)} + - (p'),, 



(32) 



(33) 



CO 

7' 



7 7 

For example, we consider the initial state 

P(0) = ^(|0) (0| + |1)(1| + |0) (1| + |1) (0|). 
The time behavior of the expectation values of position and momentum is 

(x)^. = ^^e"^^^"'^'^*cosa;t 

This is classically the motion of a damping oscillator. The time behavior of variance is 







— e' 


-(7- 


-l')t 


{p\ = 




— e~ 


-(7- 


-y)t 






— e' 


-(7- 








— e~ 


-(7- 


-7')« 



'')*) {Ap^,r + ^e-^^-^'^\l + cos' cot). 



VII. CONSTRUCTION OF LADDER OPERATORS 

In this section, we construct the ladder operators of the Liouvillian. Let L+, L_, 
be the superoperators which have the action 



L+p = p, L_p — ap, R+p — pa, R^p — 
These superoperators satisfy the following algebraic structure. 

[L_, L+] = = 1, others = 0. 
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pa 



t 



The Liouvillian can be written in terms of these superoperators. 

L = -iuj{L+L_ - R+R_) + -fL^R_ + -f'L+R+ - ^^^(L+L_ + R+R_) - 7'. 
The commutation relation between the LiouviUian and the superoperators is 



[£, L_] ^ [luj + L_ - iR^ 



We notice that {C, L^, R_) and {C, L_, R^) form closed algebras respectively, hence we 

set 

X+ = L+-R_, X_ = iR+ - 7L_ (35) 
y+ = R+- L_, y_ = 7'L+ - 7i?_. (36) 

The commutation relation between the Liouvillian and these new superoperators is 

[£, X+] = (-i(7 - 0') -tuj]x+ (37a) 



2 



[AAr_] = |i(7-7') + ^^}'^- 



(37b) 



[c,y+]^\~h-i) + ico}y+ (37c) 



2 



[A3^-] = {^(7-7')-^^}:^^- 



(37d) 



= :^+] = -(7 - 7'), others = 0. (37e) 

Therefore A'+ raises the dissipation quantum number j and raise the phase rotation quantum 
number k, y+ raises j and lowers k, lowers j and lowers k, and y- lowers j and raises k. 
There are two independent pairs of radder operators (A'+, (J^_,_, yJ). Now we determine 
the ground state. The ground state |0))) = po satisfies A". |0))) = |0))) = i.e. 

7'aVo - IPoCL^ = 0, 7'poa - lapa = 0. 
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We want to find the solution of this equation in Ti(J^). Hence we expand po on the basis 
{|m) {n\}m,n- Then this equation reduces to the equation of the expansion coefficient. We 
find 

|0))) = Ce-'^^'^^^ C = const. (38) 
Furthermore we obtain all eigenstate of the Liouvillian by acting X+,y+ on |0))). 

\m,n))) = X'lyilO))), m,neZ>o 

(39) 

£ |m, n)}} = < (7 - 7 ) - iuj{m - n) '> \m, n}}) = iim,n \m, n))) . 

Comparing the eigenvalues and the eigenstates to the previous result 

IJ'm,n = >^m+n,m-m \m , fl))) (X \m + fl , m - fl)) (40) 

Next we construct dual vectors. The dual of the ground state (((0| satisfies (((0|A'_|_ — 
(((0|3^_)_ = 0, thus (((0| oc 1. Imposing the normalization condition (((0|0))) = 1, 

(((0| = ^ - ^ V- (41) 



C V 7 
Therefore the dual vector of \m,n))) is 

^1 = , „ / , ^ (((0| A"™:^!!. (42) 

Now we have obtained the biorthogonal eigenbasis. We can write down the spectral resolu- 
tion of C and the solution of the Lindblad equation. 

00 

^ |m,n))) (((m,n| (43) 

m,n=0 

p(t) = e'^-'"*— — -— ( 1 - ^ 1 Tr {A""3^"p(0) j A'r3^?e-^'^"^" (44) 

mTI^o m!n!(7' - 7)-+" V 7^ ^ " "^^^^ + + 

This result is equivalent to the previous one. 
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Appendix A: Summary of trace class operators and Hilbert-Schmidt operators 



In this appendix we summarize the facts about the trace class operators and the Hilbert- 
Schmidt operatorai^ii^ used in this work. Let ^ be a separable Hilbert space ,^{J^) be 
the set of the bounded operators on J^, which is a C* algebra with the operator norm 
\\A\\, and {|7^)}n be a C.O.N.S of J^. The set of the trace class operators Ti(J^) = {A G 
^{J^)\ Y,n (^1 Vv4ty4 \n) < 00} is a C* algebra with the trace norm ||A||i = Trv^AU = 
^„ {n\ a/ A'^ a |n), and a two-sided *-ideal of ^{Jf). The Hilbert-Schmidt space %2i'^) = 
{A G ^{J^)\A'^A G Ti(J^)} is a Hilbert space with the Hilbert-Schmidt inner product 
{{A\B)) = TtA^B, which is well-defined because A,B e %2{^) ^ AS e %i{J^),A^ G 
12(^), a C* algebra with the norm ||v4||2 = VTrAtA, and a two-sided *-ideal of ^(Jf ). 
The relations between these three operator spaces and their norms are, 

Ti(jr) c T2(^) C ) (Al) 
\\A\\ < IIAII2, A G T2(^), ||5||2 < B G li(^). (A2) 

From (jA2]) the topology of li(^) is finer than that of T2(^). A C.O.N.S of 12(^) is 
{|m) (n|}m,nez>o- 1"^) (^1 is trace class, then Span({|m) {n\]rn,n) C Ti(^) and Ti(^) is 
dense in T2(J^). 

Appendix B: Biorthogonal system and rigged Hilbert space 

First we explain the theory of the biorthogonal basis^. In general, for a non-Hermitian 
operator A on a Hilbert space the eigenvalues of A and those of A^ are related by the 
complex conjugate. Let \mii) G be the eigenvector corresponding to the m-th eigenvalue 
Am of A, and \'mi) G be that of A'^ . Then 

A |mif) = \m \mR) , \mL) = X*ra \mL) ■ (Bl) 

But generally 7^ Imi). If the eigenvalues are non-degenerate, m ^ n ^ {^L\nR) = 0. 
If the eigenvectors are normalized, {mL\nfj) = 6m,n- We call the pair {\Tnji) , (m/^l}^ which 
satisfies {mL\nR) = 6m,n a biorthogonal system. If {\mR)}m is a complete system, the 
completeness is expressed as J2m \^r) = 1- 

Next we explain the concept of the rigged Hilbert spaced. Let Jif be a Hilbert space and 
$ be a subset of Jif which satisfies the following properties. 
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1) $ has a finer topology than that of M' . 

2) $ is complete with respect to the topology described at 1). 

3) $ is dense in M' with respect to the topology of . 

Then, the continuous anti-dual space of denoted by is an extension of the Hilbert 
space. The triplet 

$ C C (B2) 

is called the the rigged Hilbert space or Gel'fand triplet. For example, let choose I? space 
as , the Schwartz space ^ as Then is the space of tempered distributions 5^' . 

^ C C (B3) 

Now we explain the idea of the generalized eigenvector. Let $ C M' C be a rigged 
Hilbert space, A be a operator defined on and |^) e (77I e The adjoint operator 
is the operator defined on whose action is defined as 

{A^r]\i) -.^ {ri\A\i) . (B4) 

If there exist (?7a| G and A G C such that 

(AV|=A(77a|, (B5) 

that is, for any |^) e $ 

(ry,|A|0 = A(r7,|0, (B6) 

then (?7a| is called a generalized eigenvector of A corresponding to the eigenvalue of A. The 
system {(r^rnljm is called complete when 

{^ra\i) = for all m ^ 10 = 0. (B7) 

We apply the idea of the biorthogonal system to the rigged Hilbert space. Let A be an 
operator defined on {\m)}m be the system of eigenvectors of A, that is A \m) = A^ \m), 
and {(?^|}n be the system of generalized eigenvectors of A, that is, {A''n\ = A„ If the 
eigenvalues are non-degenerate, m 7^ n ^ (m|n) = 0. If the eigenvectors arc normalized, 
(m|n) — 5m,n- We call the pair {\m) ,{m\}m which satisfies {m\n) — Sm,n a biorthogonal 
system. If {|m)},„ is a complete system, {(m|}^ is also complete in the generalized meaning, 
and the completeness is expressed as |m) (m| = 1. 
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If {\m) , {m\}m is a biorthogonal system and {|m)}m is the set of eigenvectors of A and a 
complete system, {{m\}m are the generahzed eigenvectors of A. This is because any |,^) G $ 
can be expanded as |^) = Cm\fn) , Cm & C, and 

{A^n\0 = {n\ A\0 = Y^ ^-^^ ("1"^) = ^-^^ = ^« (^1^ • (^8) 

m 

Appendix C: Proof of the summation formula 

In this appendix we prove the formula introduced in section lllll This proof is due to 
Noguchi, Noumi, Yoshimoto. For p,q,r & Z>0) P ^ Q ^ f, 

[p + s - qy.{q - sy.sl 
First, changing the variable s — )■ g — s, the formula (ICip becomes 

^ (p-s)!(g-s)!s! 
We can derive the recurrence formula 

I{p, q,r - 1) = {p + q - r + l)I{p, g, r) + I{p - 1, g - 1, r - 1), (C3) 

and show that 

I{.PAA) = f^g.o, P>q- (C4) 
First we prove ( 1C2I1 in the case of r > 0. If /(p, g, g — s) = for some s < q, 

Hp, q,q-{s + l)) = {p + s + l)I{p, g, g - s) + J(p - 1, g - 1, (g - 1) - s) = 0. (C5) 

By mathematical induction, I{p, g, r) = for p > g > r > 0. In the case of r = 0, 

J(p, g, 0) = (p + g)/(p, g, 1) + /(p - 1, g - 1, 0) 

= /(p-l,g-l,0) = ■■■ = J(p-g,0,0) = 1. (C6) 
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